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All-optical active plasmonic devices are of fundamental importance for designing efficient nanopho- 
tonic circuits. We theoretically propose and numerically investigate an active plasmonic device made 
up of a nonlinear e-near-zero metamaterial slab of thickness smaller than 100 nanometers lying on 
a linear e-near-zero metamaterial substrate. We predict that, in free-space coupling configuration, 
the device, operating at low-intensity, would display plasmon mediated hysteresis behavior since the 
phase difference between the reflected and the incident optical waves turns out to be multi-valued 
and dependent on the history of the excitation process. Such an hysteresis behavior would allow 
to regard the proposed device as a compact memory unit whose state is accessible by measuring 
either the mentioned phase difference or the power, which is multi-valued as well, carried by the 
nonlinear plasmon wave. Since multiple plasmon powers comprise both positive and negative values, 
the device would also operate as a switch of the plasmon power direction at each jump along an 
hysteresis loop. 

PACS numbers: 73.20.Mf, 42.79.Vb, 81.05.Zx 



I. INTRODUCTION 

Confinement and steering of optical radiation at 
nanometer scale are fundamental issues of modern 
nanophotonics and a dominant role is played by sur- 
face plasmon polaritons 0, 0] (SPP) which are collective 
electron density oscillations coupled to external electro- 
magnetic waves Q existing at metallo-dielectric inter- 
faces. The surface charge oscillations merge the dielectric 
evanescent waves with the metal exponentially decaying 
fields so that SSPs are surface modes tightly confined at 
the interface, a fundamental ingredient for the miniatur- 
ization of optical components down to the sub-100-nm- 
size regime ]4^. Aimed at designing photonic circuits, 
controlling SPPs flow is an important topic and various 
waveguiding schemes have been proposed [H, @ together 
with a plenty of passive plasmonic devices encompas sing 
mirrors Q, lenses [H, interferometers ||, splitters [Tof 
and resonators llf. Although passive plasmonic devices 
can provide some of the essential components for sig- 
nal processing, an efficient photonic circuitry amenable 
to be integrated with nanoscale electronics evidently re- 
quires externally driven (active) plasmonic devices (l2j 
and plasmonic emitters [TH, [Hj], detectors HU and 
plasmon lasers [l?], EH nave been considered. In order to 
manipulate the plasmon flow along the path between the 
source and the detector, the dielectric refractive index 
near the metal surface has to be externally controlled 
and this can be done through an electro-optic medium 
to obtain electrically-driven plasmonic modulators fl9| . 
Electro-optic methodology is particularly intriguing since 
both optical and electrical signals can in principle be 
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processed and interconnected on the same circuit board. 
Another viable way to actively affect plasmon propaga- 
tion consists in employing a nonlinear dielectric medium 
with a large Kerr nonlinearity and this has lead to the 
observation of plasmon mediated optical bistability [2(| 
and all-optical modulation of surface plasmon polaritonic 
crystal transmission [21]. The main drawback to the use 
of the Kerr nonlinearity is that the nonlinear response is 
generally weak and, even if SSPs are accompanied by a 
large field enhancement, the required optical intensities 
are very large. For this reason, even though there is a 
considerable research effort in nonlinear plasmonics (22T - 
I2H ] , an efficient nonlinear active plasmonic device has not 
been devised yet. 

The research interest in surface plasmons has been 
further increased by the availability of metamaterials 
which are artificially structured materials whose under- 
lying subwavelength constituents are suitably tailored to 
allow the engineering of the overall medium electromag- 
netic response. As an example, it has been recently sug- 
gested that the flow of SPPs on metallo-dielectric in- 
terfaces can be efficiently (and passively) controlled by 
exploiting the metamaterial based transformation optics 
technique to devise the optical parameters of the dielec- 
tric medium [26|, [27} ■ 

The investigation of metamaterial nonlinear properties 
is particularly important in that it can lead to overcom- 
ing one of the fundamental limit of nonlinear optics, the 
fact that most of the optical materials have a relatively 
weak nonlinear response. The main idea is that the local 
electromagnetic fields of the inclusions in the metamate- 
rial can be much larger than the average value of the field 
thus producing an enhancement of the nonlinear response 
28] . Following a very different route, it has recently been 
proposed that an extremely marked nonlinear behavior 
can be observed, instead of by enhancing the nonlinear 
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response, by substantially reducing the metamaterial lin- 
ear dielectric permittivity [29| and this strategy has al- 
lowed to predict self power-flow reversing of an optical 
beam [3(| and transmission optical bistability of a slab 
of sub- wavelength thickness at very low optical intensities 
[jOLI!^- It is worth noting that, as opposed to standard 
media where nonlinearity requires macroscopic propaga- 
tion lengths (greater than a centimeter) to produce an 
all-optical manipulation, e-near-zero nonlinear metama- 
terials dramatically affect li ght just at their free-space 
interface (as discussed in Ref.[3l|) and therefore they are 
candidates to play a pivotal role in designing nanopho- 
tonic nonlinear active devices, one of which is discussed 
in the present paper. 

In this paper we theoretically investigate a nanoscale 
all-optical active plasmonic device which, operating at 
low optical intensities and being integrable in a photonic 
circuit, could be used both as plasmonic memory unit 
and as a device able to switch the direction of the power 
carried by a surface plasmon wave. The system basi- 
cally consists of a few tens of nanometers thick slab of a 
nonlinear Kerr metamaterial lying on a linear metama- 
terial substrate, the nonlinear and linear medium hav- 
ing very small dielectric permittivities of opposite signs. 
Since the substrate permittivity is much smaller than 
one, the critical angle to achieve phase-matching be- 
tween free-space waves and surface plasmons is rather 
small (of the order of few degrees) so that total reflec- 
tion of a free-space plane wave impinging onto the non- 
linear slab surface is easily observed and it is accompa- 
nied by the excitation of a nonlinear plasmonic wave at 
the slab-substrate interface. The small value of the slab 
permittivity allows the nonlinearity to produce a rather 
complex hysteresis behavior (in analogy of the transmis- 
sivity hysteresis described in Ref.[3l|) which is easily de- 
tectable by measuring the phase difference between the 
reflected and the incident optical wave. Therefore the 
system can be regarded as a nanoscale memory unit on 
a board whose state is easily read out through a reflec- 
tion experiment. Note that optical bistability from sur- 
face plasmon excitation has been investigated employ- 
ing metal-coated prisms in the attenuated total reflec- 
tion (ATR) or Kretschma nn g eometry in the presence of 
a nonlinear dielectric [33M35j ] but these setups are not 
easily scaled down to be integrated in a photonic cir- 
cuit due to the macroscopic prism dimension, a difficulty 
completely avoided by the device discussed in this pa- 
per. The nonlinear surface plasmon carries optical power 
along the slab-substrate interface and, due to the cho- 
sen extreme nonlinear regime, the overall plasmon power 
can be both parallel and anti-antiparallel to the plasmon 
wave-vector, an exotic situation generally not occurring 
in the presence of linear materials with positive magnetic 
permeability. The total plasmon power follows the hys- 
teresis slab behavior and therefore the power direction 
can be switched simply by changing either the input op- 
tical intensity or the incidence angle to force the system 
to jump to a different hysteresis state. 




FIG. 1: (a) Geometry of the plasmonic device (slab and sub- 
strate) and of the free-space incident (i) and reflected (r) plane 
waves, (b) Sketch of the electromagnetic device configuration. 
The incident and reflected waves are reported with red wavy 
lines while the arrows denoted with Si and S r represent their 
Poynting vectors. The nonlinear plasmon excited at the slab- 
substrate interface is reported with a red line and the large 
arrow denoted with J represents the overall power carried 
along the interface by the nonlinear surface mode. 



II. NONLINEAR PLASMONIC DEVICE 

We start by considering the plasmonic device consist- 
ing in a slab of thickness L = 84 nm lying on a substrate 
as reported in Figure la. Here we consider monochro- 
matic electromagnetic radiation whose free-space wave- 
length is A = 810 nm. The slab is a nonlinear Kerr meta- 
material with a very small linear dielectric permittivity 
and a viable way to produce this artificial medium is 
to repeat along the y-axis nonlinear metal and dielectric 
layers of thickness d\ and dn , respectively. If both thick- 
nesses are much smaller than A, radiation experiences a 
uniform electromagnetic response obtained by averaging 
both the linear permittivities t\ and €2 and the nonlin- 
ear Kerr coefficients xi and X2 [2^|. Therefore, since 
metal layers have negative dielectric permittivity, it is 
possible to choose d\ and di in such a way that the slab 
average dielectric permittivity is e s ; = (e,) = —0.01 (e- 

near-zero metaterial P,!^) andx (3) = (xf^ > 0. Note 
that we are considering a real permittivity e s \ and this 
can be achieved by using dielectric layers with gain to 
balance the metal losses [38[. The substrate is a lin- 
ear metamaterial with dielectric permittivity e su = 0.01 
and it can be manufactured either by means of the same 
technique employed for the slab with linear layers or by 
exploiting the method discussed in Ref. [H, H(| , consist- 
ing in dispersing silver coated spherical semiconductor 
quantum dots (with optical gain) in a homogeneous host 
material whose permittivity can be externally tuned. As 
reported in Figure la, an incident plane wave (i) is made 
to impinge onto the slab-vacuum interface (the plane 
x = 0) with incidence angle 9 and with the complex elec- 
tric field Ei belonging to the plane of incidence (parallel 
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FIG. 2: Plots of the normalized power per unit length J/ Jo 

(here Jo = -J eo/^o / i/^koX^)) anc ^ phase difference * along 
four possible ways of varying the excitation parameters, (a) 
Q = 0.006 (fixed) and 9 cr < 6 < tt/2 (variable); (b) Q = 0.015 
(fixed) and 6 cr < 9 < tt/2 (variable); (c) < Q < 0.02 
(variable) and 9 = tt/6 (fixed); (d) < Q < 0.02 (variable) 
and 9 — n/3 (fixed). In each panel, the curves plotted by blue 
empty circles represent the linear counterparts of J or ty. The 
surface Q x l obtained by Eq.((3]) is plotted in panels a3, b3, c3 
and d3 together with the critical points curve reported with 
a dashed yellow line. 

to the plane xz) which is parallel to the slab layers. A 
nonlinear plasmon mode localized at the slab-substrate 
interface (the plane x = L) can be excited only if the 
incident plane wave is phase-matched with the substrate 
evanescent waves and this happens if 9 is greater than the 
critical angle 9 cr = arcsiny 7 ?^ ~ 5.7 degrees. Therefore, 
for 9 > 9 cr , total reflection occurs since both e s i and e su 
are real and consequently |Ej| = |E r | = Qj \/%( 3 ) (where 
Q is the dimensionless electric field magnitude), the re- 
flection solely yielding a phase difference VP between the 
reflected and incident waves given by 

e** = — % (1) 
E- • p 



(see Appendix A). Evidently, the energy flow is purely 
along the z-direction since the Poynting vector is S = 
^Re (E x H*) = S z (x)e z everywhere and hence, the to- 
tal power per unit length carried by the nonlinear wave 
is 

/>oo 

J= / dxS z (x). (2) 
Jo 

We have numerically evaluated nonlinear plasmon 
modes (see Appendix A) for a number of excitation states 
(9, Q) characterized by the incidence angle 9 and the in- 
cident plane wave amplitude Q and, for each mode, we 
have evaluated the phase difference ^ from Eq.([T]) and 
the power per unit length J from Eq.([2]). In Figure 2 we 
report the evaluated J and * along four possible ways 
(denoted with a, b, c and d) of varying the parameters 
(9, Q). In situations (a) and (b) the amplitude is hold 
fixed with Q = 0.006 and Q — 0.015, respectively and 
the incidence angle is varied in the range 9 cr < 9 < tt/2 
whereas in situations (c) and (d) the angle is hold fixed 
with 9 = tt/6 and 9 = 7r/3, respectively and the am- 
plitude is varied in the range < Q < 0.02. For com- 
parison purposes, in each panel of Figure 2 reporting J 
or * we have also plotted (blue empty circles curves) 
the corresponding linear quantities i.e. J or as evalu- 
ated for a linear slab of dielectric permittivity e s i. Note 
that, in each situation reported in Figure 2, both J and 
^ are multi-valued function of the excitation parame- 
ters or, more precisely, this quantities generally assume 
more than one value for a given excitation state (0, Q). 
This implies that the considered system displays multi- 
stability controllable through the incidence angle and the 
input optical intensity, and detectable by measuring J or 
"F . From Figure 2, it is evident that the values of J and 
^ lie on at most three different branches (red, green and 
black lines) and, most importantly, that the branches 
presents points at which they abruptly stop (see, e.g., 
state F of panels al and a2, states M and K of panels 
bl and b2, state B of panels cl and c2 and states F and 
H of panels dl and d2). The fact that different branches 
abruptly stop and that they are not mutually continu- 
ously connected imply that the multistability discussed 
in this paper is supported by a mechanism fundamentally 
different from that (feedback effect) yielding the standard 
optical bistability. 



III. PLASMONIC MULTISTABILITY 

In order to discuss the physical mechanism supporting 
the obtained multi-stable phenomenology it is convenient 
to consider the electromagnetic matching conditions at 
the interface x = L which yield 



As schematically sketched in Figure lb, total reflec- 
tion produces a nonlinear localized mode at the inter- 
face x = L which is a transverse magnetic (TM) field 
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where Q x l — Q X (L) and Q z l — Q Z (L) are the values 
of the x and z dimensionless electric field components 
at the interface x = L on the slab side (see Appendix 
A). For 9 and Q z l assigned, Eq.Q provides the allowed 
values of Q x l (electric field component normal to the 
slab-substrate interface). Note that for a linear slab the 
first term of Eq.© reduces to e s iQxh this implying that 
only one value of Q x l is allowed. In the four panels a3, 
b3, c3 and d3 of Figure 2 we plot the surface Q x l ob- 
tained from Eq.Q in the range Q x l > (without loss 
of generality) and we note that such surface is folded. 
This implies that, due to the effect of the nonlinearity, 
the electromagnetic boundary conditions generally allow 
more than one electromagnetic configuration at a given 9 
and plasmon amplitude Q z l, the various configurations 
having different values of Q x l- Each point of the surface 
corresponds to a nonlinear plasmon mode and therefore, 
in panels a3, b3, c3 and d3 of Figure 2, we have also re- 
ported the various state curves of plasmon modes used 
to evaluate J and "J of the corresponding panels in the 
Figure, using the colors to establish the correspondence. 
The red branches of both J and \t of Figure 2 practi- 
cally coincide with their linear counterparts (the empty 
circles curves) apart from the region close to the breaking 
points so that we conclude the green branches to be of 
pure nonlinear origin and arising from the surface fold- 
ing which is therefore the crucial element producing the 
multistability of J and <I>. It is worth noting that folding 
takes place if the surface admits critical points where the 
partial derivative of Eq.© with respect Q x l vanishes, 

1. e. 

est + \ (9QIl + QIl) = 0, (4) 

which represents the curve of critical points reported with 
a dashed yellow line in panels a3, b3, c3 and d3 of Figure 

2. From Eq.Q we readily deduce that critical points can 
appear only if the field amplitudes \Q x l\ and \Q z l\ are 
comparable with y/—e s i and this physically means that 
linear and nonlinear contributions to the overall electro- 
magnetic slab response must be comparable. This is an 
extreme nonlinear condition [29j made here possible by 
the fact that the slab permittivity e s ; is much smaller 
than one. It is worth noting that, in a conventional 
nonlinear medium where the permittivity is of the or- 
der of unity, the above folding would required an op- 
tical intensity so large to effectively prevent the whole 
discussed phenomenology to take place. We therefore 
conclude that the multistability discussed in this paper 
can uniquely be observed by means of e-near-zero non- 
linear metamaterials. To get further physical insight, 
note that the constitutive relations yield D x (L,0) = 
e [e.j + \ (3Q 2 xL + Q 2 zL )] Q xL /^W\ so that Eq.flJ) co- 
incides with the condition dD x (L,0)/dQ X L = 0. There- 
fore the plasmon states below and above the critical 
curve are characterized by dD x (L,0)/dQ X L > and 
dD x (L, 0)/dQ X L < 0, respectively and therefore char- 
acterized by very different electromagnetic features (the 



curve of critical point playing the role of a separatrix 
between the two electromagnetic regimes). This is par- 
ticularly evident from panel d3 of Figure 2 where states 
H and F are just below and above the critical curve, 
respectively (and hence their electric amplitudes are al- 
most equal at x = L) and nonetheless their correspond- 
ing values of J and 'J of panels dl and d2, respectively, 
are very different. The existence of two distinct electro- 
magnetic regimes explains both the different values of J 
and "J along the red and green branches (bistability) and 
branch breaking since, from Figure 2, it is evident that 
the points at which the branches of J and abruptly 
stop correspond to a plasmon mode whose Q x l and Q z l 
lie just either below or above the yellow critical points 
curve (state F of panels al, a2 and a3, state M and K 
of panels bl, b2 and b3, etc.) 



IV. DEVICE MEMORY AND POWER 
SWITCHING FUNCTIONALITIES 

By combining the discussed multistability with the 
branch breaking it is easy to prove that the considered 
system displays a complex hysteresis behavior allowing 
to regard the system as a device with structured mem- 
ory functionality. Suppose, for example, that the incident 
plane wave is such that 9 = tt/3 and Q = 0.001 so that, 
from panels dl, d2 and d3 of Figure 2, it is evident that 
a state of the red curve is excited since for Q < 0.004 the 
green branch has no corresponding allowed states. By 
increasing the input intensity (and hence Q) the state G 
is reached and surpassed toward the state H. Once H 
is reached, if the input intensity is further increased the 
system has to undergo a sudden jump to the state / be- 
longing to the green branch since the red branch stops 
at H. From the state I the input intensity can now be 
decreased toward the state F and therefore along such a 
backward path the attained values of J and are dif- 
ferent from the forward path, i.e. hysteresis occurs. Ev- 
idently, by decreasing the input intensity from the state 
F the system undergoes another jump to reach the state 
G, thus coming back to the red branch and closing the 
hysteresis loop. The system hysteresis phenomenology is 
however far more rich since the incidence angle can also 
be varied. Suppose, for example, that the incident plane 
wave is such that 9 — tt/6 and Q = 0.015 so that, from 
panels bl, b2 and b3 of Figure 2, it is evident that the 
state D belonging to the green branch is excited since 
the red branch has no corresponding allowed state. If 
the angle is increased, by reasoning in analogy with the 
just discussed hysteresis loop, the system can be brought 
to the state K (passing through the state N) and hence 
forced to jump to the state L since K is a breaking point 
of the green branch. From L the system can be driven 
to M by decreasing the incidence angle and eventually it 
can be forced to jump down to the state N thus closing 
the "angular" hysteresis loop. 

The number of possible hysteresis loops is infinite since 
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FIG. 3: (a) Global multi-valued surface of the normalized 

power per unit length J/ Jo (here Jo = U eo//W(2fcox' 3 '))- 
(b) Global multi- valued surface of the phase difference 9 . 
Dashed lines represent the branches of J and 9 considered 
in Figure 2 and analogously for the states reported with cap- 
ital letters. 



9 and Q can independently be varied and, in order to 
provide a global description of the system multistable 
behavior we report in panels (a) and (b) of Figure 3 the 
values of J and <F, respectively, for all the possible 9 and 
Q in the ranges 9 cr < 9 < n/2 and < Q < 0.02. In 
Figure 3 we also plot, through dashed lines and corre- 
sponding colors, all the branches of J and \1/ considered 
in Figure 2 together with all the states denoted with cap- 
ital letters. The surfaces representing J and 'F exhibit 
more than one sheet, each one admitting breaking lines, 
and the mutual position of the just discussed hysteresis 
loops G-H-I-F and N-K-L-M is particularly evident. If 
we choose to operate uniquely with the red and green 
sheets, it is evident that the system is suitable to record 
a binary information which can be easily read since a 
single measure of >F (or J) allows to know whether a red 
or a green state is excited. Moreover, the state of the 
memory unit can be simply altered by driving the sys- 
tem through a breaking line thus forcing a switch of its 
operating sheet. Note also that the system memory func- 
tionality can even be suitably tailored to satisfy external 
requirements (i.e. specific optical intensities or angles 
imposed by a possible nearby circuit environment) since, 
for each plasmonic state, a hysteresis loop can be found 



which starts and ends at the considered state. Suppose, 
as an example, that an hysteresis loop containing the red 
state A of Figures 2 and 3 is required. Once the state 
A is excited, one can increase the field amplitude Q to 
reach the state B (which is a red sheet breaking point) 
and to force the system to jump to the green state C. 
From panel cl and c2 of Figure 2 or from Figure 3 it is 
apparent that by solely varying the field amplitude the 
system can never be brought back to the state A. How- 
ever, one can drive the system to reach the state E, fix 
the field amplitude Q and increase the angle 9 to reach 
the state F (see panel al and a2 of Figure 2 or Figure 3) 
which is a breaking point of the green sheet. Therefore a 
jump can be induced to make the system reach the red 
state G and, after reducing the angle, to come back to 
the initial state A thus closing the hysteresis loop. 

In addition to the above discussed multistability and 
related memory functionality, we note that, from Figures 
2 and 3, both positive and negative signs of J are allowed. 
Specifically, the red branches of J are always positive 
whereas the green branches can be both positive and neg- 
ative. Bearing in mind that the red branches (where they 
exist) almost coincides with the slab linear behavior, we 
conclude that the extreme nonlinear regime also intro- 
duces the possibility of an overall plasmon power which 
is antiparallel to the wave- vector. In fact, since the slab 
and substrate permittivities e s i and e su are negative and 
positive, respectively, the Poynting vector is locally an- 
tiparallel and parallel to the z-direction (which coincide 
with the plasmon mode wave- vector) and hence the total 
power can be both positive and negative. For the red 
branches of J of Figure 2 (or 3) the positive contribution 
to J coming from the evanescent field in the substrate 
overcomes the slab negative contribution thus produc- 
ing an overall positive J. However, the plasmon modes 
belonging to the green branches are generally character- 
ized by a larger amount of electromagnetic energy stored 
within the slab since their amplitudes Q x are generally 
larger than their red counterparts as it is evident from 
panels a3, b3, c3 and d3 of Figure 2. Therefore the neg- 
ative contribution to J of the green modes can overcome 
the positive evanescent contribution thus yield an over- 
all negative plasmon power. As an example, states A 
and E of panels al of Figure 2 are excited through the 
same input intensity and incidence angle and, nonethe- 
less, J(A) > and J(E) < 0. Obviously, for small values 
of 9, the extension of the evanescent field thorough the 
substrate is so large to always yield a positive value of J 
(see the positive portion of the green sheet of J in Figure 
3a). 

The fact that the overall plasmon power, for a given ex- 
citation state (9, Q), can be both parallel or anti-parallel 
to the wave vector also allows to regard the considered 
system as a device able to switch the direction of the 
electromagnetic energy flowing along the slab-substrate 
interface. In fact, once a state is excited with a given sign 
of J, the energy flow direction switch is easily performed 
by externally changing the incident plane wave and driv- 
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ing the system through a sheet breaking line whose cross- 
ing forces the system to jump to the other sheet where 
the sign of J is generally different from the initial one. 



V. DEVICE FEASIBILITY AND DESIGN 

The feasibility of the above described active plasmonic 
device is evidently related to the actual availability of 
both slab and substrate metamaterials exhibiting the 
above mentioned features. In order to prove that such 
media can be actually synthesized, we have chosen d\ = 2 
nm, d 2 = 5 run, L = 84 nm, ex = -28.81 + 9.78i, 
e 2 = 11.51 - 3.911 Xx 3) = 3.16 x 1CT 16 m 2 /V 2 and 

(3) 

X 2 = 0. The parameters of medium 1 coincide with 
those of silver 14 1|. characterized by a very large nonlin- 
ear susceptibility [42J , with the imaginary part of the per- 
mittivity corrected by the layer size effect (since d\ — 2 
nm) (43|. Medium 2 is a linear dielectric with gain 
(to compensate the metal losses) and its permittivity €2 
can be achieved, for example, by embedding, in a poly- 
mer background, a bi-dimensional array of colloidal II- VI 
quantum dots which are properly electrically or optically 
pumped [44| through ultra-violet light. The value of e 2 
here considered is obtained through the Maxwell Gar- 
nett mixing rule where = 10.2 and cqd = 10 — 15« 
[44j are the background polymer and quantum dots per- 
mittivities, respectively and / = 0.29 is the quantum 
dots volume filling fraction; besides both the polymer 
and the quantum dots have nonlinear susceptibility much 
smaller than xi so that the requirement | X2 1 ^ Ixil is 
assured. Note that, even though the actual fabrication 
of the considered ultra-small structure is hindered by the 
very small thickness/ width (i.e. di/L) ratio of the metal- 
lic grooves, the required array of metallic nanowires can 
be manufactured through the general method of Super- 
lattice Nanowire Pattern Transfer (SNAP) (H-HJ which 
is a novel deposition technique particularly suitable for 
producing ultrahigh-density arrays of aligned metal and 
semiconductor nanowires and nanowire circuits. Once 
the array of silver nanowires is deposited onto the sub- 
strate, the vacancies between the metallic grooves can be 
filled with the required quantum dots dispersed in a suit- 
able liquid monomer (by using, for example, the standard 
spin-deposition process) which can be eventually poly- 
merized through UV curing. The effective (average) per- 
mittivity and nonlinear susceptibility of the considered 
sample are e s i — {d\e\ + d 2 e 2 ) / {d\ + d 2 ) = — 0.01 + 0.00H 
andx (3) = (d lX ( i ) +d 2 X2 ) )/(di+d 2 ) = 9 x lir 17 m 2 /V 2 . 
As far as the substrate, we choose here to exploit the just 
proposed scheme of dispersing silver (for which €A g — 
—28.81 + 9. 78*) coated spherical semiconductor quantum 
dots (the same colloidal II- VI quantum dots as in the slab 
for which cqd = 10 — 15z) in a homogeneous host ma- 
terial (for which eh = 10.2). Exploiting the electrostatic 
treatment of Ref. and choosing p = 0.19630 for the 
fraction of the total particle volume occupied by the inner 
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FIG. 4: (a): Comparison between the phase difference 9 eval- 
uated through full wave simulations (blue circles) for 9 = n/6 
and < Q < 0.025 and the corresponding phase difference 9 
of panel (c2) of Figure 2 (red and green solid lines), (ul), (u2), 
(u3) and (vl), (v2), (v3): profiles of the plasmonic squared di- 
mensionless electric field components as a function of x (slab: 
< x < 84 nm; substrate: x > 84 nm) corresponding to 
the states U and V reported in panel (a). The electric field 
profiles have been evaluated at y ~ di, e.g. close to the 
metallic-dielectric interface. 



quantum dot core we obtain for the effective permittiv- 
ity of the overall inclusion = —17.87 + 0.05i. Ex- 
ploiting the standard Maxwell Garnett mixing rule and 
choosing / = 0.04625 for the volume filling fraction of 
the inclusions dispersed in the hosting matrix, we obtain 
for the effective dielectric permittivity of the substrate 
e s „ = 0.01244 + 0.000011 



VI. FULL- WAVE SIMULATIONS 

In order to test the behavior of the above discussed 
active plasmon device we have performed 3D full-wave 
simulations by means of a finite-element solver (see Ap- 
pendix B) in which we have used the numerical values of 
the parameters just chosen in the above structure design. 
In panel (a) of Figure 4 we report the comparison be- 
tween the phase difference W evaluated through full wave 
simulations (blue circles) for 9 = 7r/6 and < Q < 0.025 
and the corresponding phase difference \& of panel (c2) 
of Figure 2 (red and green solid lines). We note that 
full- wave simulations confirm the existence of the above 
discussed multistability and branch breaking. For val- 
ues of Q far away from the breaking value Q = 0.01, 
full-wave simulations reproduces the predictions of the 
homogenized approach. On the other hand, discrepan- 
cies appear for 0.01 < Q < 0.015 and, specifically, the 
breaking point of the upper branch predicted by the full- 
wave simulations is at Q — 0.015 instead of Q = 0.01 of 
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the homogenized theory. The origin of the discrepancies 
and of the breaking point shift is due to the existence of 
a further plasmonic resonance arising from the fact that 
a transverse magnetic plane wave (polarized in the xz 
plane, see Figure la) is made to impinge onto a metallo- 
dielectric medium layered along the y-direction. There- 
fore an additional y component of the electromagnetic 
field of plasmonic origin arises thus departing the layered 
medium behavior from that of its homogeneous counter- 
part. In order to support this reasoning we plot in panels 
(ul), (u2) and (u3) and panels (vl), (v2) and (v3) of Fig- 
ure 4, the profiles of the squared dimensionless electric 
field components as a function of x corresponding to the 
states U and V reported in panel (a), respectively. Such 
squared electric field profiles have been evaluated near 
the metallic-dielectric interface (y ~ d\) where the addi- 
tional plasmonic states are mainly localized. The state 
V does not display significant discrepancies between full- 
wave simulations and homogenized approach prediction 
(see panel (a) of Figure 4) since, as it is evident from 
panels (vl), (v2) and (v3) of Figure 4, both the x and 
z electric field components are much larger than the y 
component. On the other hand, the discrepancy accom- 
panying the state U are due to the fact that the y electric 
field component of panel (u2) is comparable with the z 
component of panel (u3). 

The intensity of the incident plane wave is / = 
(1/2) v^o/MoQVx' 3 ^ which, for the amplitude range Q < 
0.02 (where multistability occurs for e s ; — —0.01) and for 
the above effective nonlinear susceptibility x > yields 
/ < 1.17 MW/cm 2 , which are intensities smaller than 
those normally required for observing the standard op- 
tical bistability. However it is evident that the smaller 
\e s i\ the smaller the intensity required for observing the 
device multi-stable behavior so that very smaller inten- 
sity thresholds (of the order of W/cm 2 ) are very likely to 
be attained. 



VII. CONCLUSIONS 

In conclusion we have theoretically considered an all- 
optical active plasmonic device characterized by a com- 
plex multi-stable behavior allowing the device, in prin- 
ciple, to be used both as a memory unit or as a switch 
of the plasmonic power direction. Since the transverse 
dimension of the proposed device active region (slab) is 
smaller than 100 nm and the system is controlled through 
free-space coupling, the device could be easily integrable 
on a photonic circuit where it can be possibly intercon- 
nected with other nanophotonic and/or microelectronic 
devices. As opposed to standard situations where the 
nonlinearity mediated plasmonic flow steering is ham- 
pered by the weakness of the nonlinear response and 
the short propagation distances, the nanophotonic de- 
vice proposed in this paper fruitfully exploits the nonlin- 
ear response since the linear dielectric permittivities are 
very small and, consequently, nonlinearity rules the elec- 



tromagnetic behavior just at the slab-substrate interface. 
We believe that our methodology offers a new way for ex- 
ploiting optical nonlinearity to actively control plasmonic 
propagation and hence for the designing of future active 
plasmonic devices with complex functionalities. 



Appendix A: Nonlinear Plasmonic Modes 

We have numerically investigated the excitation of 
nonlinear plasmonic modes at the slab-substrate in- 
terface by solving the slab nonlinear Maxwell equa- 
tions together with the electromagnetic matching con- 
ditions at x = and x = L. Specifically, 
the TM electromagnetic field within the slab is 



E 



[Qx{x)e x + iC 



^i(ko sin )z 



H 



^JJT Q [Q y {x)e y ]e^^ z / ^) (where k Q = 2tt/X 
and Q x , Q y and Q z are, without loss of generality, di- 
mensionless amplitudes) and it satisfies Maxwell equa- 
tions with the Kerr nonlinear constitutive relation D = 
e {e.iE + x (3) [(E • E*) E + § (E • E) E*] } for < x < 
L, or 



1 dQ z 
ko dx 

1 dQ y 
ko dx 

Q y sin 9 



sin 9 — Q y 



1 



tsl + \ {3Ql + Ql 



Qr, 



(Al) 



The free-space (x < 0) incident and reflected plane 
waves are E s = E s (e y x k s /fco) exp (ik s • r), H s = 
E s \J e /fi e y exp (ik s • r) (where s = i,r, Ei, E r 
are the amplitudes of the two waves and k^ = 
ko (cos8e x + sin#e 2 ), k r = fc (— cos8e x + sm6e z ) are 
the wave vectors) whereas the evanescent field within the 
substrate (x > L) is 



E = E„[ sin 9e x — iy sin 6 — e su e 



-k$-\J sin 2 9 — e su (x — L) „i(ko sin 0)z 



H — E v 



3 — feo y sin 2 — e su (x — L) ^i(ko sin 9)z 

(A2) 



where E p is, without loss of generality, the real field am- 
plitude and 9 > 9 cr is assumed. At the slab input (x = 0) 
and output (x — L) interfaces we require the continuity 
of the tangential electric field component (E z ) and of 
the normal displacement field component (D x ) (the con- 
tinuity of the tangential magnetic field component H y 
evidently follows from the third of Eqs. ljAlj )). In order 
to evaluate a nonlinear plasmonic mode we specify the 
amplitude E p so that the matching conditions at x = L 
yield 



Qz 



sin 9 — e. 
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+ Q 



zL) 



QxL = 



, sine 



(A3) 



where Q x l = Q X {L) and Q z l = Q Z (L) are the values of 
the x and z dimensionless electric field components at the 
interface x = L on the slab side. After Q x l and Q z l are 
obtained from Eqs. (IA3l ). Eqs. (|Al[) are numerically solved 
all the way to the input face at z = where the match- 
ing conditions allow to evaluate the amplitudes of the 
incident and reflected waves, Ei and E r , which turn out 
to have the same absolute value. Note that, by combin- 
ing the first and the second of Eqs. (|A3[) . Eq.Q is easily 
obtained. 



Appendix B: Full- wave simulations 

The phase difference VP, plotted as blue circles in Fig- 
ure 4(a), together with the profiles of the plasmonic 
squared dimensionless electric components, reported in 
Figure 4, are obtained from a 3D full-wave simulations 
using the COMSOL RF Module. In order to perform 
the simulations corresponding to the setup of Figure 1, 
we have considered an integration domain containing a 



unit cell consisting of a vacuum layer (in the region 
— 160 < x < nm, < y < 7 nm, < z < 1 nm), 
a metal layer characterized by e% — —28.81 + 9.78i, 
Xi 3) = 3.16 • 10- W m 2 /V 2 (in the region < x < 84 nm, 
< y < d\ = 2 nm, < z < 1 nm), a gain medium with 

(-2 = 11.51 — 3.9D, %2 = ( m the region < x < 84nm, 
di < y < d = 7 nm, < z < 1 nm) and a homoge- 
neous substrate with e su = 0.01 (84 < x < 1704 nm, 
< y < 7 rim, < z < 1 nm). Periodic boundary con- 
ditions have been imposed at the faces perpendicular to 
the y— axis (at y — and y — 7 nm) for achieving an 
infinite perfect tiling of the system whereas periodic Flo- 
quct boundary conditions have been imposed at the faces 
perpendicular to the z— axis (at z = and z = 1 nm) for 
adapting the simulation domain to the z field periodicity. 
Scattering boundary conditions have been imposed at the 
faces perpendicular to the x— axis (at x = —160 nm and 
x = 1704 nm) for both exciting the incident plane wave 
and to absorb the waves scattered by the structure. The 
simulation was performed by using a rectangular nonuni- 
form mesh adapted to the structure design and the con- 
sidered radiation excitation. 
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